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Probably the most fundamental processes for biological systems is their ability to create themselves 
through the use of cell division and cell differentiation. In this work a simple physical model is 
proposed for biological cell division. The model consists of a positive ionic gradient across the cell 
membrane, and concentration of charge at the nodes of the spindle and on the chromosomes. A 
simple calculation, based on Coulomb's Law, shows that under such circumstances a chromosome 
will tend to break up to its constituent chromatids and that the chromatids will be separated by 
a distance that is an order of thirty percent of the distance between the spindle nodes. Further 
repulsion between the nodes will tend to stretch the cell and eventually break the cell membrane 
between the separated chromatids, leading to cell division. The importance of this work is in 
continuing the understanding of the electromagnetic basis of cell division and providing it with an 
analytical model. A central implication of this and other studies is to give theoretical support to 
the notion that cell division can be manipulated by electromagnetic means. Requirements on the 
ingredients of more sophisticated models for biological cell division will also discussed. Copyright 
2006-2010 by the author. 
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INTRODUCTION 

Understanding mitosis and cytokinesis, the pro- 
cesses by which biological cells divide, has pro- 
found implications to practically all areas of 
medicine and biology. It also exposes the very 
important physical question: how does a com- 
plex, self organized, animate structure develop 
from mechanistic-inanimate constituents. An un- 
derstanding of cell division and its manipulation 
could lead, e.g., to the ability to inhibit the growth 
of tumors, the ability to restart growth in defective 
or damaged cells, as well as many other exciting 
prospects. Unfortunately, though there has been 
significant progress in understanding cell division's 
mechanical parts, [ll-Q a clear picture of the major 
underlying physical mechanisms involved does not 
exist. Such an understanding is the most efficient 
method to harness the potential in manipulating 
cell division. 

In this report a simple mechanistic approach is 
presented for cell division. In this approach it is as- 
sumed that a charge gradient is created across the 
cell membrane through protein interactions and 
that positive ions such as Ca"*"^ and Na^ enter 
through these ion channels into the cell's spindle 
structure, a network of protein based tubes under- 
lying the innards of the cell during cell division. 
Though this is a feasible construct it is not 



claimed here to be the actual mechanism of cell 
division but rather as a test for the required elec- 
tromagnetic interactions. Out of the four funda- 
mental forces known in nature it is only the elec- 
tromagnetic force that can play a major role in 
small scale biology in general and cell division in 
particular. Charge gradients play a significant role 
in other biological systems such as nerves. 



Theory 

Let us discuss the potential energy of the sys- 
tem. To allow the model to be transparent, only 
the bare minimum cell components are included: 
a cell membrane, two spindle centers with charge 
qi = qi, and two chromatids with charge (72 =93- 
To a given charge qi, with i = 1,2,3,4, a one di- 
mensional position Xi{t) and velocity Vi{t) are as- 
sociated. The coordinate system is chosen such 
that, xi = — X4 and X2 = —2:3 initially. A diagram 
of this configuration is given in Fig. [TJ The ini- 
tial velocity is taken to be zero for all four charges. 
These assumptions will insure that xi = —X4 and 
X2 — —X3 for all times. The rational for choos- 
ing to work on a one-dimensional problem goes 
beyond the attempt to present the most trans- 
parent model. The structure of the dividing cell 
would define a preferred spatial direction along the 
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line joining the two spindle centers, hence a one- 
dimensional model is appropriate for the position 
of the four charges. 




Cell Membrane 



FIG. 1: Components of the minimal cell. A cell 
membrane, two spindle centers, qi and 54, and two 
chromatids, q2 and g3, which carry the genetic infor- 
mation. Figure not to scale. 

The classical electrostatic contribution to the en- 
ergy can be expressed as 

i^.. = ^4 + 1 + 8,,,.^), (1) 

with I — Xi ~ xi, d = X3 — X2, and D is the inverse 
dielectric constant. Here it will be assumed that 
the cell can be modeled as water and Z) = 80 is 
taken. The symmetry of the system about the ori- 
gin was used to arrive at Eq. ([T]) . The fact that the 
electrostatic potential energy is positive relative to 
the charges at infinity comes about from the repul- 
sive interaction between the positive charges, since 
the system would prefer to be in a lower energy or 
a more separated state. 

A well known result for the electrostatic field 
inside biological media and other neutral solution 
is the Debye length. This calculation shows that 
for a neutral solution with free positive and neg- 
ative ions distributed uniformly, the electrostatic 
field due to an ion drops off exponentially once 
the field it generates is in the vicinity of opposite 
charges. The reason the notion of a Debye length 
does not invalidate the model proposed here is that 
the charge distribution in the dividing cell is not 
neutral, nor does the mitotic spindle allow the ho- 
mogeneous distribution of ions throughout the cell. 
It has been known for several decades that charge 
pumping into the cell and mitotic spindle accom- 
panies cell division, (see for example Ref. Q p. 



263). Hence, The charge density in the dividing 
cell is not properly accounted for under the as- 
sumption of a neutral solution. One should also 
remember that the existence of positive charges in 
the cell's media does not reduce the electrostatic 
repulsion force between charges qi and in Fig. [1] 
it increases it. 

At this point it is instructive to find the minima 
of the electrostatic energy, as this will illuminate 
on how the charged chromatids separate the ge- 
netic code of the cell. If we assume for simplicity 
that Iq = 10 /im and constant, the minimum of 
Eq. ([Ij is given by the implicit condition for do, 

l/dg = 16do/o/('o ~ ^0)^1 91 = 12- This con- 
dition can be solved to give do = 3.6148 /im, i.e., 
if we start with a configuration of chromatids that 
is close together they will tend to separate to a 
distance that is of the order of magnitude of the 
initial cell size. Of course based on Eq. ([T|) I itself 
would tend to increase as well, which will finalize 
the separation of the cell. However, as we shall 
soon see changes in I occur on a much larger time 
scale than changes in d. 

The next ingredient in this presentation is the 
deformation of the cell membrane. This is a cru- 
cial step both for the initial division of the cell, and 
equally important, in the capability of the daugh- 
ter cells to separate from each other. Here the ap- 
proach is to assume that the cell membrane is a lin- 
ear catastrophic spring, i.e., that it resists stretch- 
ing by a force that is proportional to its pertur- 
bation in length and that given some catastrophic 
condition the spring will break, and the force be- 
tween its two ends will be zero, representing the 
separation of the two cell membranes. This force 
is applied to the spindle centers, qi, q^ but not to 
the chromatids, 52, qs, which are assumed to be 
small in comparison to the cell size, free of the 
tensions of the cell membrane, and approximately 
not affected by the nuclear membrane. The poten- 
tial energy associated with a linear spring can be 
written as 

Est = ^Hxi{t) ~ ^M? + ~ ^M?. 

(2) 

with fc the spring constant. 

The catastrophic condition is taken to be a 
stretching of the cell membrane to twice its length. 
When this happens it is assumed that £'57- = 0. 
Under the assumptions presented so far, maximum 
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stress will be experienced by the cell niembrane at 
the position a; = 0. Hence it will be assumed that 
the cell membrane will break at x = 0. The ac- 
tual separation of the two daughter cells involves 
joining of opposing cell membranes around x = 0. 

The movement of the chromatids inside the cell 
is important to the process of separation. The de- 
tails of the separation of a chromosome into its 
two chromatids probably involves the breaking of 
molecular bonds, for example by charge substitu- 
tion. However for the rest it will be assumed that 
the details of separation of the chromatids can be 
modeled as a separate process and will be ignored. 

The details of the various forces the chromatids 
experience in the process of division are again rela- 
tively complex. The drag force acting on the chro- 
matids, i.e., on q2 and ^3, is 

Fd,^ = -iv,{t) (3) 

with i = 2, 3, 7 is a dissipative parameter. As shall 
be shown shortly, this parameter is important to 
the behavior of the chromatids. 

RESULTS 

We can get a relatively quick estimate for the 
number of charges involved in the process by con- 
sidering a simpler two charge problem. Let us de- 
fine the average force, F, from the initial combined 
stage, Xi, and the final separated stage, x/, 

/ Fdx= Fdx^-{EEs{xf)-EEs{x^)), 

J Xi J Xi 

(4) 

with F = —dEss/dx, including for the moment 
only the electrostatic interaction, and x is chosen 
as the position of one of the charges. Now, we 
can associate with this average force an average 
acceleration, a, obeying 

Xf -x^ = ^d{tf - U) (5) 

with ti and tf the initial and final times. 

The change in the electrostatic potential energy 
for two identical charges with charge q is given by 
q^{l/xf — 1/xi). Combining with Eqs. (jH) and ^ 
it is found that 

9 2(x f ~ x.i)mXiX f „ , , 



with TO the accelerated mass, m can be estimated 
as the mass of a sphere of water lO/im in radius, 
TO = 10~^f;. If we use Xi — 5/iTO, Xf — 10/iTO, 
and tf — ti = s ds a, typical time scale for the 
separation phase, we find from Eq. a charge 
that is equivalent to about 10 Ca+^ atoms. This 
order of magnitude estimate should be viewed as 
a feasibility test for the proposed theory. 

With the aide of a numerical simulation we can 
much better calculate the dynamics of the sys- 
tem. For simplicity it will be assumed that due 
to the structure of the spindle the two halves of 
the splitting cell can be represented by the mass 
of the dominant charges toi and TO4. The sys- 
tem is started with four charges at rest, situated 
at —xi{t = 0) = X4^(t = 0) = 5/im, —X2{t = 0) = 
x^{t = 0) = 10 nm, with charge qi = = lOe, 
12 ~ (l3 — 2e, where e is the charge of the pro- 
ton. The main accelerated mass is composed of 
the charged ions but mostly of other cell materi- 
als contained in half the cell, toi — m4 = 10~^ g. 
The mass of the chromatids is estimated as the 
mass of water occupying a tube of length ten fim 
and radius ten nm, TO2 — TO3 = 10~^^ g. The 
parameter k in Eq. Q is estimated by the arbi- 
trary requirement that at the catastrophic condi- 
tion the force due to the surface tension would 
equal the electrostatic force between qi and (74. 
Hence, k = 10"^° g/s^. 

Using these values and 7 = 1~^° g/s for the 
parameters determining the system, the dynamic 
evolution is simulated. The system is started from 
the initial condition and an average acceleration is 
calculated, then it is evolved by a At, new positions 
and velocities are calculated using the calculated 
accelerations, etc. It was found that a time step 
of 1 ns produced convergence to within graphical 
resolution, i.e., reducing At further did not change 
the graphical appearance of the solution. 

Results from such a calculation are presented in 
Fig. [2] and Fig. |3l Since the aim of this calcula- 
tion is to present a single mitotic event it was as- 
sumed that after the catastrophic event, when the 
cell has elongated to twice its original length, both 
Est = and Ees — 0. The time when this condi- 
tion is imposed is indicated by an arrow in Fig. [5] 
The differentiation of the original cell is clearly ev- 
ident in Fig. [5J Note that the original distance 
between the two chromatids is not visible on the 



4 




10 
5 

-5 
-10 
-15 
-20 







3 4 
t [s] 



FIG. 2: Position as a function of time for the 
system described in the text. —xi{t = 0) = X4{t — 
0) — 5 jj.m, ~X2{t = 0) — xz{t — = 10 nm, gi — 
q4 — lOe, q2 = ~ 2e, with e the charge of the 
proton, mi — m4 — 10~* g, m2 — ma — 10~^^ g. 
k = 10"* g/s^ 7 = 1"^'' g/s. Both the cases of At = 
1 ns and At = 10 ns are graphed and produce identical 
results to within graphical resolution. 



scale of the figure, while the eventual separation of 
the chromatids is of the order of magnitude of the 
original cell size. 
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evident from the figure that due to the relative 
closeness of the chromatids and their relatively low 
mass, their response times are much shorter than 
the those of the relatively massive poles. Hence, 
the dynamics in the initial stages will be dominated 
by that of the chromatids. This relates back to the 
previous conjecture that the chromatids act as ini- 
tiators of cell division. Another important aspect 
to be gathered from Fig. [3] is that the parameter 
7 plays a role in the movement of the chromatids 
and hence the dynamics of cell division. For low 7 
the chromatids will undergo more oscillations than 
for high 7 where the chromatids slowly move to 
their separated configuration. The case presented 
in Fig. [3] is mostly steady relaxation though there 
is some overshoot in the separation between the 
chromatids. 

Finally, it is interesting to investigate the rela- 
tionship between the time of separation and the 
number of charges involved. The condition for sep- 
aration that will be used here is that the distance 
between qi and 54 is double its original distance. 
The time that it takes the cell to reach this point 
is denoted as t2 ■ In Fig. |4l ^2 is given as a func- 
tion of the charge of qi and 94 which is denoted by 
q. It should be observed that the for high q the 
electrostatic effect from the charges (72 and (73 be- 
come negligible and i2 cx 1/g as expected for two 
charges. For low q the effects of the charges q2 and 
^3 become more important and ^2 deviates from 
l/q dependency. 
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FIG. 3: Blowup of Fig. [2l Used to show the earlier 
response of the chromatids. The drag parameter 7 will 
determine the type of relaxation the chromatids will 
experience. Both the cases of At = 1 ns and At = 10 ns 
are graphed and produce identical results to within 
graphical resolution. 

The initial details of the separation are presented 
in Fig. [21 which is a blowup from Fig. [21 It is 



CONCLUSIONS 

In summary, a transparent model was presented 
exhibiting cell division involving charge separation. 
It was shown that through the resulting electro- 
static interaction the cell will separate into two 
symmetric parts, each containing one pole and one 
chromatid. 

Probably most important is that cell division is 
inherently an electromagnetic phenomena. Funda- 
mentally there is no other choice. Hence, it should 
be possible to manipulate cell division by electro- 
magnetic fields. For example, in the treatment of 
cancer it may be possible to inhibit tumor growth 
by creating electromagnetic fields that disfavor cell 
division. The converse is also true, it may be pos- 
sible to stimulate tissue growth by electromagnetic 
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FIG. 4: Time of separation. t2, as a function of 
the charge of gi and §4 which is denoted by q. For 
high q values the charges 92 and 53 become negligible 
and t2 <xl/q as expected for two electrostatic charges. 
For low q values the effects of the charge of q2 and §3 
become more important and t2 deviates from the 1/q 
dependence. 

fields that favor cell division. 

The presentation given here took the most ba- 
sic approach, so as to be most transparent. Topics 
such as the regulation of charge content through- 
out the cell, which were mostly ignored here, are 
critical to the further understanding of and forth- 
coming observations about cell division. 



After the formulation of theses ideas a related 
paper was found. Q In it the effects of high- 
frequency alternating electric fields on cancer tis- 
sue growth are studied. For frequencies in the 
range 100 — 300 kHz it was found that alternat- 
ing electric fields have an inhibitory effect on the 
ability of the cancer cells to divide. Other evi- 
dence suggests that low frequency alternating elec- 
tric fields with frequencies less than approximately 
1 kHz can promote tissue growth in some cases. [6| 
While not directly supporting the work presented 
here these studies strongly suggest the manipula- 
tion of cell division is possible with electromagnetic 
fields and this in turn confirms that cell division is 
inherently an electromagnetic phenomena. 

The author would like to acknowledge the con- 
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Electronic address: lfaraggi@ph ysics.utexas.edu ] 
[1] J. S. Hyams and B. R. Brinkley, ed., MITOSIS: 
Molecules and Mechanisms (Academic Press Inc., 
1989). 

[2] M. Glotzer, Curr. Opin. Cell. Biol. 9, 815 (1997). 

[3] L. J. Gagliardi, Phys. Rev. E 66, 011901 (2002). 

[4] P. K. Hepler, in MITOSIS: Molecules and Mecha- 
nisms, edited by J. S. Hyams and B. R. Brinkley 
(Academic Press Inc., 1989), chap. 7, pp. 241-271. 

[5] E. D. K. et al.. Cancer Research 64, 3288 (2004). 

[6] C. A. Basset, Clin. Plast. Surg. 12, 259 (1985). 



